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We discuss a simple, semiclassical scattering theory for spin-dependent 
transport in a many-terminal formulation, with special attention to the four 
terminal device of Johnson referred to as spin transistor. 
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1. Introduction 

Transport in ferromagnets and ferromagnetic multilayers has has become 
a popular subject of research since the discovery of the giant magnetoresistance 
(GMR) effect, but has been studied long before. Tedrow and Meservey [[jj showed 
that the tunneling current in superconductor/insulator/ferromagnet junctions is 
spin dependent. Under an applied bias, the current of spin-up and spin-down elec- 
trons injected from a ferromagnet into a normal metal is different. This causes a 
non-equilibrium magnetization or " spin accumulation" which gives rise to an addi- 
tional boundary resistance, the "spin-coupled interface resistance" fl [I. Johnson 
and Silsbee M detected this spin accumulation by two ferromagnetic contacts to 
a normal metal. Spin injection can be the physical basis of new devices, like the 
so-called spin transistor 0, p5J. A "pedagogical model" of this device is shown 
in Fig. [11(a). It consists of a normal metal film sandwiched between two ferro- 
magnets. The magnetizations of the two ferromagnets are either aligned parallel 
or antiparallel. A spin-polarized current is injected from the first ferromagnet 
Fl to a drain contact P connected to the normal metal. The second ferromag- 
net F2 is attached to a floating voltage probe. The polarized current creates 
a non-equilibrium magnetization in the normal metal, which splits the chemical 
potentials of the spin-up and spin-down electrons. Just as Fl can be seen as a 
spin polarizer, F2 can be seen as a spin detector, at which the chemical poten- 
tial aligns with that of the spin-up (spin-down) electrons in the normal metal 
in the parallel (antiparallel) configuration. The voltage measured at F2 is used 
to determine the impedance difference between the parallel and the antiparallel 
configuration, which is related to the spin-coupled resistance. The experiments 
are actually carried out with an additional normal metal counter electrode as in 
Fig. ^(b). Spin-dependent transport is affected by the properties of the normal 
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Fig. 1. Schematic picture of the spin transistor for (a) the pedagogical model and (b) 
the actual configuration. 

metal. Spin-flip scattering in the normal metal mixes the two spin channels and 
thus decreases the spin accumulation. The spin-orbit interaction || appears to 
limit the performance of the spin transistor. We discuss in this paper the multi- 
terminal Landauer-Biittiker scattering theory of transport for a magnetic system. 
This formalism provides an alternative view on the physics of the device |Q. |8| 
which in the diffusive regime is equivalent to previous two-terminal theories by 
Johnson and Fert and Lee ||], but we demonstrate that it can be extended to 
other regimes as well. For technical details we refer to pi . 

2. The spin transistor in the Landauer-Biittiker formalism 

In a two-terminal configuration the current from a contact A to a contact B 
is / = G(/.ia — A*s)e, where A/i = /ia — Hb is the difference in chemical potential 
between the two contacts. The linear response conductance G is related to the 
transmission probability by the Landauer formula |7]]: 
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where t m(jna i is the transmission amplitude of an electron from mode m with spin 
a in lead A to mode n and spin a' in lead B, all states being at the Fermi energy. 
Here we assume for simplicity that the magnetizations are collinear, i.e. parallel 
or antiparallel to each other. A network theory of magnetic multi-terminal devices 
with general magnetization is presented in [nTl . 

Buttiker extended the two-terminal Landauer formula to a many terminal 
device by summing over all contacts H. If Ta^b is the total transmission prob- 



ability from contact A to B: 
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T^ B = ]T^X> mff ,„ ff f, (2) 
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where N a is the number of modes in contact A and Ng the number of modes in 
contact B. The total current in contact A is given by 
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Referring to Fig. EYa), we call the transmission probability from Fl to F2 
T\. Those from the ferromagnetic contacts Fl and F2 to the drain contact of the 
paramagnetic film are taken to be the same T 2 . When the net current into F2 is 
zero, the potential of F2 can be obtained using Kirchhoff's law: 

V S = V-^—. (4) 

V s depends on the choice of the potential zero, here taken to be here in the drain 
contact reservoir. On the other hand, Z s , the impedance between Fl and P, only 
depends on the potential difference V. A constant current / is driven from Fl 
into the drain contact of the paramagnetic film, which according to Eq. (0) is 
equal to 

e 2 
I = --T 2 {V + V S ). (5) 

The impedance Z a between Fl and P is therefore 
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For the (pedagogical) system with two fully polarized ferromagnets without spin- 
flip processes the transmission possibility between the two ferromagnets is zero 
when they are aligned antiparallel. The impedance Zf p is then simply the recip- 
rocal conductance between Fl and P 

7 ap _ h 1 
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The difference between the parallel and antiparallel configuration is 
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where T p is the transmission probability to F2 in the parallel configuration. 

In the presence of spin-flip scattering the transmission probability T x no 
longer vanishes. The difference in impedance between the parallel and antiparallel 
configuration now becomes 

rpP _ rpAP 
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Johnson B implicitly assumed that the current into the drain contact did not 
affect the impedance difference. This can be modelled by an infinite resistance to 
the drain or T2 = 0: 
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Eq. ([!(]) is equivalent to Johnson's result, which diverges in the absence of spin- 
flip scattering T 1 AP — » 0. 

Instead of the ill-defined V s , Johnson measured a voltage V s with respect 
to a normal counter electrode N, c.f. Fig. [j](b). The transmission probabilities 
from Fl and F2 to the normal counterelectrode N are called Tfin and Tp2N, 
respectively and T/vp is the transmission probability from N to P. Again setting 
the net current into F2 and N to zero: 
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Assuming that the resistance to the drain contact is relatively large, we can con- 
sider two limits. When the normal counterelectrode N is positioned far from the 
drain contact P both transmission probabilities to the drain contact, T 2 and T^p, 
vanish and the impedance of the system reduces to 
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The impedance difference is 
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Comparing Eq. (R4h with Eq. (1Q), we expect the four terminal configuration to 
be equivalent to the three terminal configuration when Tp 2 N *C Tfin- However, 
the impedance difference has a maximum for Tp 2 ^ = 0: 
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which is increased by a factor 4/3 compared to the pedagogical model. 

When the normal counterelectrode N is positioned close to the drain contact 
P, Tnp is no longer negligible and Tpuy and Tp 2 N are small. The impedance 
difference 
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now depends on the ratio of Tpuv and Tp2N- It has a maximum for Tpix — Tp2N 
for which it reduces to Eq. (|10|). Apparently, the impedance difference is affected 
by the choice of position of the normal counterelectrode N . 

A microscopic calculation of the various transmission probabilities involves 
bulk, spin-dependent interface and spin-flip scattering processes. We present here 
a first attempt based on a semiclassical calculation of a weakly scattering slice of 
a disordered material. 

3. Sheet impurity scattering with spin-flip 

In the effective mass approximation the single-electron states at the Fermi 
energy Ep are described by the Schrodinger equation in a normal metal 

n 2 



2m* 



V z + V(r) -Ep)I+ H S9 (f) + H so (r) 



^ T (f) ^ _ 



0. 



(17) 



where m* is the effective mass. We consider here the effect of a thin sheet of short- 
range impurities with a scalar scattering potential V(r) = ^2 a ^a5{z)5{p — p a ), 
where p a gives the transverse position of the scattering center, and j a gives 
the strength of the scatterer. Spin-flip scattering can be induced by spin-orbit 
scatterers Q 

HS ° = 4^2 - ' [(V V(r)) x (-ih V)] , (18) 

where ~a r is the vector containing the Pauli spin matrices, or interaction with 
magnetic impurities with fixed random spin direction (spin glass): 
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Here Sa is the spin of the paramagnetic impurity and J ex is the local exchange 
integral. 

We wish to compute the Green function |13, ft% H 
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where ?R t t — {pr,l, zr,l) and zr > 0, Zl < are located in the right and left 
leads. Without perturbations 
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The matrix of transmission coefficients is related to the Green's function by (l3| 
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For isolated interfaces the following optical theorem holds jlj] : 
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We are interested in the transport properties averaged over impurity configura- 
tions. The averaged Green function reads 
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The conductance can be calculated in two ways. By using Eq. (|23|) and the 
relation (E4) between the Green's function and the transmission coefficients: 
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Alternatively, the conductance can be calculated diagrammatically via the trans- 
mission probabilities, which can be written in general as 
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Diffuse scattering and spin diffusion are described by the second term of Eq. (|2€ 
which is governed by the (reducible) vertex correction W . 

These formally exact relations can be evaluated for the Born approxima- 
tion which is valid for low impurity densities nm = Nm/A and weak scattering 
strength. With an average strength of the scatterers 7 = J2a^a,a/NjB, and a 
mean square value j 2 = ^2 a 7 2 a /Nm, the self-energy reads in the Born approx- 
imation 

ft 2 
S B = n IB j~i—k F (r] + 7i s f), (27) 
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where r\ is the spin-conserving scattering parameter 
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'7s/ = 'ylf+^s? is the spin- flip scattering parameter caused by spin-orbit scatterers 
with u so — h 2 kp/Am 2 c 2 

2n m f m*u sol V 

and by paramagnetic impurities 
n iR frn*J x 
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Assuming ReS B < ImS B and inserting Eq. (Elj) into Eq. (25) 
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To the first order in 77 the conductance is given by G/Gq = 1 — 2(77 + ?/ s /) , where 
Go = (2e 2 /h)(Akp/4:ir) is the Sharvin conductance. In order to calculate the 
spin diffusion, we need the transmission probability for spin-flip scattering. To 
the lowest order in the scattering parameters the transmission probability matrix 
in spin space is obtained by summing (\tt p\ 2 ) over the incoming and outgoing 

states. The spin-flip probability T s f under transmission is 

T„ - ^\n„- (32) 

whereas the spin-conserving transmission probability T is 

4. Finite thickness 

4-1- Relevant length scales 

Next to the geometrical parameters of a sample, i.e. the sample cross sec- 
tions A and length L, several characteristic length scales govern the transport 
properties. In the (quasi-) ballistic regime the conductance is dominated by the 
contact resistance. In the diffuse regime transport is limited by scattering at bulk 
impurities. An obvious parameter is the impurity scattering mean free path £. The 
spin-flip mean free path t a j is the average length an electron travels before it flips 
its spin. Including scattering at spin-orbit scatterers or paramagnetic impurities, 
does not only decrease £ s f — vpT s f, where T 8 f is the spin-flip scattering time and 
vf the Fermi velocity, but also reduces the mean free path £~ x — Vp 1 ^^ 1 +t~}). 
In the diffuse regime the relevant parameter is not £ s f, but the spin diffusion 
length l s f — y/££ s f/6 over which the spin-accumulation persists |I5| . 

4.2. Two-terminal conductance 

The transport properties of samples with finite thickness can be readily ob- 
tained from the results of the previous chapter in two different limits. When the 
mean-free path is much larger than the current path, the system is in the quasibal- 
listic limit and the results above can be carried over directly. In the quasi-ballistic 
regime, L <C £, we can make the connection of the scattering mean free paths and 
the microscopic parameters N(n + rj s f) = L/2£ and Nn s f = L/2£ s f. The con- 
ductance is G qb (L)/Go = 1 — 2L/£. The spin-flip probability for transmission 
through a thin slice of bulk material follows from Eq. (p2|): 



8 



and for transmission without spin-flip 

In the diffuse regime, where L > f or JV(i) + r) s f) ^> 1, the conductance is 
G d f(L)/Go = 4£/3L. The spin-resolved two-terminal transmission probabilities 
in the diffuse regime, on the other hand, can be obtained from the scattering 
properties of a thin slice by following Schep et al. p| . It is shown in fill that 
this procedure is equivalent to solving the diffusion equation |l5|, |) with the 
following general solutions for the spin-average ~p = Az + B and the spin-splitting 
Sfi of the chemical potentials: 

S/j,(z) = Cexp(—z/l s f) + Dexp(z/l s f). (36) 

The integration constants A, B, C, D have to be determined by the boundary con- 
ditions. The transport properties are easily obtained from the chemical potentials. 

5. Spin transistor 

We can now collect the different results to obtain expressions for the spin 
transistor. We approximate the transmission probabilities in the many terminal 
configuration by those in the two-terminal configuration. This is not unproblem- 
atic because we do not know the current paths. Complication are for example 
spreading resistances at contacts, which reflect the fact that electrons do not al- 
ways travel straight between the two contacts. However, in a thin film geometry 
like the spin transistor, we expect the corrections to be not too large for a qual- 
itative study. It should be kept in mind that the geometric parameters, i.e. the 
cross section A and the current path length L, may somewhat deviate from the 
geometrical measures of the sample. 

5.1. Quasi-ballistic regime 

Let us consider first the quasiballistic regime and assume that the contacts 
are made from strong ferromagnets in which the minority carrier density vanishes 
and the majority density equals that of the normal metal island. In that limit the 
performance of the spin-transistor is determined by the spin-flip scattering length 
£ s f > £. From Eqs. (B4|) and (|35|), the transmission probabilities are 
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The transmission probability to the drain contact can be cither quasi-ballistic 
or Ohmic, depending on the distance Li between the ferromagnetic and drain 
contact. If the drain channel is quasi-ballistic (Z2 *C I) 

Akl f L 2 \ 

T2 = ^[ 1 -t)- (39) 



The impedance difference for the pedagogical model is to lowest order in L\ 
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This equation also holds when transport to the drain is diffusive. 

5.2. Diffuse regime 

Let us now consider tunnel contacts with a spin-injection efficiency which 
can be expressed in terms of G^ , the tunneling conductance for spin a, respec- 
tively the density of states at the Fermi energy N a (E F ): 
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We will now consider a system consisting of a normal metal of thickness Ln sand- 
wiched between two ferromagnetic tunneling junctions, over which a voltage Afi/e 
is applied. Defining the average tunneling conductivity G F = 2Gf Gf /(Gf +Gf ) 
we can solve the diffusion equations, with the result for the transmission coeffi- 
cients: 
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where gN = Akpi/Sn. In the limit T 2 — ► the impedance difference becomes 
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has a maximum for G F = 0: 
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which agrees with [|9| in the tunneling limit. 

By an analysis of his experiments Johnson found l s f = 1.5±0.4/xm ||. The 
impedance of the spin transistor is approximately 25 /if2 in the parallel configura- 
tion and the measured impedance is 3 /iS! for the 1.6 /iin thick film. The largest 
impedance difference predicted for the optimum value |/?| = 1, using Eq. (46) and 
the above value found for l s f, is Z s Ad « 3 x 10~ 2 il /im 3 , which is about three 
times smaller than the measured difference. This problem might be caused by 
contamination which reduces the contact area. 
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6. Conclusions 



We analyzed the spin transistor by a semiclassical scattering theory of trans- 
port resolving some inconsistencies which arise from a two-terminal approximation 
for a many-terminal device. The performance of the spin transistor is limited by 
the spin-diffusion length in the diffuse regime and the spin-flip scattering length 
in the quasi-ballistic regime. Naturally, we many expect a better transistor action 
in clean samples. 
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